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ABSTRACT 

The advent of satellite -to-satellite tracking has provided a means for recovering 
atmospheric refractivity profiles and for retrieving certain atmospheric param- 
eters such as pressure and temperature from radio occultation data. Profiles 
calculated in this manner might be used to monitor world-wide weather patterns. 

The purpose of the present investigation is to demonstrate the feasibility of re- 
covering these parameters from one-way range rate between two earth orbiting 
spacecraft during occultation of the tracking signal by the earth's lower atmos- 
phere, One of the spacecraft, the ATS-6 {Applications Technology Satellite), is 
in a geostationary orbit in the equatorial plane. The other, NIMBUS-6, is in a 
near-circular almost polar orbit of about 1100 kilometers. 

The tracking data is inverted by an integral transformation (Abel transform) 
to obtain a vertical refractivity profile above the point of closest approach of 
the ray connecting the satellites. The only restriction concerning the refractivity 
function is that it is spherically symmetric. 


After accounting for ionospheric effects and water vapor in the lower atmosphere, 
either by independent measurement or regression techniques, pressure and 
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temperature distributions can be obtained from values, of dry refracttvity using 
the hydrostatic equation and perfect gas law. 

Two methods are investigated for recovering pressure and temperature param- 
eters. In the first method refractivity (which is proportional to density in a dry 
atmosphere) is numerically integrated to obtain pressure. Temperature is then 
calculated from pressure and dry refractivity. In the second method the param- 
eters of a linear temperature profile mth height (base temperature and lapse 
rate) over a given height region are solved for in the least squares sense using 
values of dry refractivity. Weighted apriori information (such as a standard 
atmospheric temperature profile or previous radiosonde or satellite radiometer 
measurements) can also be included. Pressure is then obtained from the tem- 
perature profile and dry refractivity at the base of the atmospheric layer. 

An error analysis is performed in order to assess the sensitivity of random and 
bias type errors in measured range rate, satellite position and velocity, and 
height upon recovered parameters. Results show that recovery is much more 
sensitive to satellite velocity errors than to satellite position errors, in par- 
ticular the position errors of the geostationary satellite ATS-6 where an error 
of 10 kilometers in each position component effects an error of about 0.2% in 
recovered refractivity. An error in the height of the geoid of 5 meters effects 
an error in recovered refractivity on the order of 0. 10%. 
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An example is given demonstrating recovery of parameters from radio, occulta- 
tion data obtained during satellite-to-satellite tracking of NIMBUS-6 by the 
ATS-6 satellite. Pressure recovery by integration of dry refractivity agreed 
more closely with radiosonde data with a maximum deviation of about 3. 5%. The 
least squares procedure for estimating temperature, with and without apriori in- 
formation, yielded a temperature profile closer to radiosonde data with a maxi- 
mum deviation of about 1.5%. 
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RECOVERY OF REFRACTIVITY PROFILES AND PRESSURE AND 


TEMPERATURE DISTRIBUTIONS IN THE LOWER ATMOSPHERE 
FROM SATELLITE-TO-SATELLITE OCCULTATION DATA 

1.0 MTRODUCTrON 

The advent of satellite-to- satellite traeMng [1] has provided a means for re- 
covering atmospheric refraetivity profiles and retrieving certain atmospheric 
parameters such as temperature and pressure from radio occultation data. 

Radio occultation of a tracking signal between two satellites occurs whenever 
the signal passes through the atmosphere (Fig. 1), When this occurs, the ray 
is bent. In the lower atmosphere the bending is opposite to that in the iono- 
sphere,* In addition, the signal is delayed as it travels through the atmos- 
phere, but advanced as it travels through the ionosphere. The first effect is 
due to the variation in the index of refraction across the signal path which re- 
siiLts from changes in atmospheric temperature, pressure, and relative humid- 
ity, The second effect is due to the fact that the index of refraction is greater 
than unity in the lower atmosphere, but.less than unity in the ionosphere (phase 
index of refraction). At 2 MHz which is the frequency of the carrier signal for 
satellite-to -satellite tracking in this analysis, the ionosphere has less of an 
effect than the lower atmosphere, especially for rays that pass close to the 
earth. Therefore, both the bending and delay along the path combine to cause 

"For this analysis the lower atmosphere is less than 100 kilometers. 
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LOW SATELLITE ORBIT^ 



NOTE: 

FIGURE IS IN PLANE CONTAINING BOTH SATELLITES AND 
THE CENTER OF THE EARTH. 



Figure 1. Occultation Geometiy 



the apparent range between the two satellites to be larger than the line of sight 
range and the range rate along the ray path to be larger than the range rate along 
the line of sight. 

The mathematics of profile inversion are well-documented [2] and planetary 
atmospheres have been studied using microwave radio occultation techniques for 
Mars ( [3] , [4], [5]) and Venus { [6] , [7], [8]). Morrison, Ungar, and 
Lusignan [9] have proposed a system of satellites (one master satellite and a 
number of slave satellites) for the purpose of measuring pressure levels in the 
atmosphere. The satellites .are placed in a common circular orbit such that the 
propagation path between the master and each of the slave satellites probe dif- 
ferent levels in the atmosphere. Kliore [10] analyzed the feasibility of such a 
system and showed how changes in the observable quantity (phase path change) 
are related to changes in such meteorological parameters as density , tempera- 
ture, pressure, and humidity . 

The basic problem in this analysis is the following. Given the effect of the 
atmosphere upon one-way range rate (or equivalently doppler) between two 
earth-orbiting satellites (one in geostationary orbit at approximately 36,000 
kilometers, the other in a low orbit of about 1100 kilometers), what is the re- 
fractivity profile which produced this effect? And, can temperature and pres- 
sure distributions be retrieved from the refractivity profile ? 
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The .answer to the first question is straightforward. With a knowledge of the 
position and velocity ofboth satellites .during radio occultation,* the measured 
range rate can be used in an integral transformation (Abel transform) to obtain 
directly a vertical reffactivity profile of the lower atmosphere above the point ' 
where the ray path comes closest to the earth [11] . The refraetivity N is as- 
sumed to he spherically symmetric and is related to the index of refraction by 

N = N(r) = [n(r) - 1] lO’® (1) 

where n(r): is the index of refraction for the lower atmosphere at the geocentric 
radius r. 


Refraetivity, however, is a function of temperature, pressure, and relative 
humidity, and is given by ( [12] , pg, 7) 

77.6 P 

N = + (3.719) 10^ 

T 



where P is pressure in millibars, T is temperature in degrees Kelvin, and e is 
the partial pressure of the water vapor in the air (in millibars) given by ( [13] , 
pgs. 343, 358) 


e 


7-5 (T - 273.15) 

— ^ (6 1 1 ) - 273 . 15 ) 

100 


(3) 


where Rj^ is relative humidity in percent. 


*The position and velocity of each satellite during occultation are determined from tracking data taken 
during nonocculted time periods using an orbit determination program. The state vectors of the satel- 
lites at a given epoch (time point close to occultation) are solved for in the program, then propagated 
through occultation to provide position and velocity of the satellites at the data times. 
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From (2) it can be seen that recovering temperature and pressure from re- 
fractivity is not as straightforward as recovering refr activity from range rate 
data. Account must be taken of ionospheric effects and water vapor which con- 
tributes to the wet component of refractivily, the second term on the right-hand 
side of (2). If these effects can be properly accounted for, either by independent 
measurement or regression techniques , all that will remain will be the dry 
component of refractivity , the first term on the right-hand side of (2). This 
quantity is proportional to density in a dry atmosphere. 

In this analysis two methods for obtaining temperature and pressure distributions 
from dry refractivily are investigated. In both cases the equation of hydrostatic 
equilibrium and perfect gas law are assumed to hold. In the first method dry re- 
fractivity is numerically integrated to obtain the pressure. Temperature is then 
obtained from pressure and dry refractivity using the first term on the right- 
hand side of (2). In the second method the parameters of a linear temperature 
profile with height in a given atmospheric layer (base temperature and lapse 
rate in that region) are obtained from values of dry refractivily using nonlinear 
least squares. Apriori information (such as a standard atmospheric temperature 
profile or previous radiosonde or satellite radiometer data) can also be included. 


*By selecting data above a certain height (for example 5 to 7 kilometers), the water vapor effects in the 
lower atmosphere can be neglected. At 2MHz and for rays passing close to the earth, the change in the 
contribution of the ionosphere to the total path length (as well as the rate of change of path length) of 
the signal should be rather small. For these reasons a simple model for ionospheric density in the 
upper atmosphere should suffice. For rays passing higher up in the atmosphere a dual frequency sys- 
tem can be used to remove ionospheric effects. 
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Pressure is then obtained from the recovered temperature profile and the dry 
refractivily at the base of the layer. 

The advantage of obtaining pressure by numerically integrating dry refractivity 
is that no assumption is made concerning the temperature structure. In a dry 
atmosphere refractivity is a direct measurement of density. In order to re- 
cover the pressure accurately, however, a sufficient number of points over a 
wide range of heights must be available for the numerical integration. If there 
is sufficient reason to believe that temperature is linear with height in a given 
region, then not as many points would be required to obtain the least squares 
solution for the temperature profile. 

There is a similarity in the mathematical techniques of recovering temperature 
from dry refractivity and recovering temperature from radiation measurements, 
as well as in the structuring of temperature as a piecewise linear function with 
height. Smith [14], [15] models the lower atmospheric temperature with two 
lapse rates. In regard to the subject of temperature retrieval, Rodgers [16 ] 
has written an excellent comprehensive survey paper. 

An advantage afforded by the geometrical configuration of the satellite orbits 
in this analysis (one in geostationary orbit and the other in a noncoplanar orbit) 
is that different pressure levels in the lower atmosphere can be probed with just 
one satellite pair. In addition, for a fixed Inclination of the nongeo stationary 
satellite orbit, occultations will occur at about the same latitude in both the 
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Figure 2. A' 


3.0 TRACKING MEASUREMENTS 


For the ATS- 6/NIMBUS- 6 Tracking, Experiment [l],the signal travels a 4-way 
path as shown in Figure 3 (from [1]). The range measurement consists of the 
.time, in .seconda.betw:e.en ,a zero, .crossing of the^ground transmitted signal (.6 GHz 
carrier phase modulated by a 100 KHz sidetone) and a zero crossing of the ground 
received signal (4 GHz carrier having two subcarriers each phase modulated by 
the 100 KHz sidetone) after the signal has travelled from ground to ATS, to 
NIMBUS, back to ATS, and then to ground [1]. The range rate measurement 
consists of the time in seconds required to count a fixed number of cycles of the 
4-way Doppler frequency plus a fixed offset frequency. 

For inversion purposes the measurement consists of the one-way range rate be- 
tween ATS-6 and NIMBUS-6 which is obtainable from the 4-way Doppler data.* 

4.0 ATMOSPHERIC EFFECT 

The atmospheric effect upon range is the difference between the electrical range 
along the ray path and the geometric strai^t line distance between the two sat- 
ellites (absence of an atmosphere). The atmospheric effect upon range rate is 

’"The 4-way Doppler output for ATS-6 Satellite-to-Satellite tracking can be approximated by [17] 

fd = [ajfi + 32 (rj + f2)] 

where = measured average Doppler frequency 
fj = uplink frequency (6137.85 MHz) 
c = speed of light 

k, a j , 32 are scalar constants determined by equipment frequency multiplication 
(k~.336,ai =0.76,32 = 1.0) 
il = average one-way range rate ATS to ground site 
T 2 = average one-way range-rate ATS to Nimbus. 
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ATS-6 



ROSMAN, N.C. 

tracking station 


Figure 3. Tracking Measurements (from [1]) 


the derivative of this difference, Referring to Figure 4, the electrical dis- 
tance, Rg, along the path NDCBA is given by 

Rg = I n(r)ds (4) 

where n(r) is the index of refraction for the lower atmosphere and ds is an 
incremental length along the ray path. 


The atmospheric effect upon range and range rate is therefore 


AR = Re 

-Rg 

(5) 

AR ~ Rg 

-Rg 

(6) 
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An example of the atmospheric effect upon range rate for an actual ATS-6/ 
NIMBUS-6 occultation pass on 15 April 1976 can be seen in Figure 5. 

i 

5.0 HERGLOTZ-WIE CHERT INVERSION TECHNIQUE 

This technique was originally used to determine velocity-depth profiles in the 

earth from the travel times of seismic waves between known positions of a 

source and a receiver located on the earth's surface. It has been adapted for 

the radio occultation problem by Phinney and Anderson [11] (see Appendix) . As 

1 

used here the technique provides a direct means for transforming satellite-to- 
satellite tracking data obtained during occultation into a profile of refractive 
index versus radius. 

The data used In the radio occultation problem is one-way range rate between 
the two satellites (obtainable from the 4-way Doppler data described in Section 
3,0). The method requires as input two different range rate residuals. One 
residual is used to determine the angle between the ray path and the radius 
vector of each satellite which is necessary to determine the impact parameter 
(Section 6.0). The other residual is used in the integral inversion formula to 
determine the refractive index at the minimum probing distance. 

The main advantage of the Herglotz-Wieehert method compared to, say, a 
model fitting approach, is that the only assumption made concerning the analytical 
form of the refractiviiy is that it is spherically symmetric. Thus, no modelling 


errors are introduced. 
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Referring to Figure 4 and following Phinney and Anderson, let r^ be the radius 
which defines the outer limit of the lower atmosphere (for r > ro the index of 
refraction is unity). A ray parameter, p, called the impact parameter of the ray, 
is defined as - 

p = n(r)rsini (7) 

where n(r) is the index of refraction at the radius r and i is the angle of incidence 
between the tangent to the ray path and the radius vector. 

By Snell's law p is a constant along a ray, and in the absence of an atmosphere, 
is the radius of closest approach of the .ray. With an atmosphere present p is 
perpendicular to the ray asymptote. We can also write 

p = To sin io = HpFp (8) 

where the subscript p refers to the turning point of the ray. 

Define a new variable 17 

r? = n(r) r = (9) 

san 1 

which equals r^ when r = ro and p when r = rp. The electrical range along the 
ray path between points B and D (which is proportional to the total phase shift at 
the carrier frequency) can then be expressed as a function of p 

(t? 2 - drj (10) 

where ds is an increment along the ray path. 


R(P) 


= 2 y* nds = 2y 
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But, the range along the ray path between B and D can be written as the sum of 
the range in the absence of an atmosphere, R^, plus the range residual due to 
an atmosphere, R^ . 

R = Ru + Ra (11) 


where R^ = 2pcotio. 


Using (7) through (11) and the Abel transform (see Appendix), the refractive 
index np at the minimum probing radius rp is given by 


Hp = exp 


'--f 


0 /p' 

cosh"-^ — 

V’? 


1\ /dR. 


dp 


dp 


( 12 ) 


Using (4) and the fact that the ray paths from B and D to the satellites are 
straight lines (ionospheric bending Is negligible), Rg^ can be written as 




= f n(r) ds - A/Ri-p2 - x/r^-p^ 


= Re - \/ Ri - - \/ Rn - P^ 


(13) 


Differentiating (13) with respect to p gives the residual dRa/dp in terms of the 
measured one-way range rate Rg between the satellites 


d]^ 

dp 




1 1 

( , )* + p 
\dp/ 

VRa-p" \/R^-p"J 


(14) 
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The integration in.;(12) is carried out numerically using Gaussian quadrature. 

In order to use. (M)- in the jntegrand^of (12), it is nec.essary.to determine p as a 
function of t. For this purpose a least squares polynomial is fit to t versus p 
and dt/dp obtained from the polynomial. 

Refractivity is obtained using Tjp in equation (1). 


6.0 DETERMINATION OF IMPACT PARAMETER 

The impact parameter p can be obtained iteratively using as input the ephemeris 
of both satellites and a range rate residual [18]. 

In Figure 4, which is-’in the plane. containing both satellites and the center of the 
earth, the angle i// is given by 


cos - 


■ ^NR 


iv, 


N‘ 


(15) 


where is the low satellite velocity vector projected onto the plane of both 
satellites and the earth center, and is the projection of along the ray 
path. 


The measured range rate Eg can be expressed as 

“ ^NR ^AR ( 16 ) 

where Var is the projection of the high altitude satellite velocity vector along 

/ 

the ray path (the dot product of the velocity vector with a unit vector along the 
ray path at the satellite and directed away from the low satellite) and is 
similarly defined for the low satellite. 
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Combining (15) and (16) we have 


^ = COS"^ 


Re - V 
_ ivnI 



(17) 


The quantity (Eg - V^j,^) in (17) is a range rate residual. 

The angle j3jg is given by 

~ ’^'e (18) 

where is the angle between Vjj and the line of sight and ajq is the angle 
between the radius vector of the low satellite and the line of sight. Both 
these angles are determined from the ephemeris of the satellites. 


The impact parameter is given by 

P Rfj sin/3jsj 


(19) 


The iteration proceeds as follows. The projection of the high satellite velocity 
vector along the line of sight is input as a first guess for in (17) and the 
angle calculated using the measured range rate and IVnI { iVj^l can be 
determined from the position and velocity of the satellites). Next is obtained 
from (18) and used in (19) to determine a value for p. Using this value for p and 
(radius vector of the high satellite) and the ephemeris, a new value for 
is obtained which is then input to (17). The process is continued until conver- 
gence is reached. 
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From the above it is seen that two range rate residuals are required as input in- 

» 

the Herglotz-Wiechert Inversion technique. The residual (R^ - ^ar) used in 

> 

(17) to obtain the impact parameter p iteratively, while the other residual dRg / 
dp, given by (14), is used in the integral transformation of (12) to obtain the 
index of refraction Up at the minimum probing radius Tp. 

7.0 IONOSPHERIC EFFECTS 

At S-band frequencies (2 GHz carrier tracking signal between ATS and NIMBUS) 
and particularly for rays passing close to the earth, -.the effects of the ionosphere 
upon range .and range. fate are small compared to the effects of- the lower' atmos- 
phere [19], The way in which the ionosphere enters info the occultation cal- 
eolations is in the computation of the path length of the signal through that region 
(100 kilometers to 1000 kilometers). For rays close to the earth's surface, the 
change in the contribution of the ionosphere to the total path length is small. 

Its effect upon range rate, therefore, should be negligible. Figure 6 shows the 
effects of the lower atmosphere and ionosphere upon range, range rate, 
and bending angle for an assumed exponential model for lower atmospheric re- 
fractivity and a Chapman model for the electron density In the Ionosphere with 
parameters corresponding to a total electron content of 10^^ electrons/meter^ . 

The exponential refractivlty .is given by • 

Nj = Nq exp(-h/H.j.) = (nj - 1)10® (20) 
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IONOSPHERIC EFFECTS 
RANGE ERROR AR, ■ IQ-"' (AR, IN METERS) 
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Figure 6. Comparison of Lower Atmospheric and Ionospheric Effects on 
Range, Range Rate, and Bending Angle 


with No = 330. 3 

H-p =6.6 kilometers (scale height) . 
iij = index of refraction, 
h = height in kilometers 


The Chapman model used in the Analysis is given by 


Ne = 


N 


max 


exp 


|(1 - z - e-z) 


with Njjjgjj = (0.3098) • 10^^ ele ctrons /meter ^ 

= (h - h^)/Hi 

h„^ = 300. 73 kilometers (height of maximum electron content) 

. ’• ^ ' M 

•Hj = 78.11 kilometers (scale height) 

‘ I ‘ ■ 

h = height in kilometers 


( 21 ) 


Since the basic measurement involved in satellite- to-satellite tracking is the 
doppler shift of the tracking signal/ithe phase index of refraction is used for ray 

} I 

tracmg through th’e ionosphere, and the refractivlty is given by [20] 


Ni = - 


40.3 N,(h)‘ 


10 ^ 


(nj - 1 ) • 1 0.^.. 


( 22 ) 


with nj = phase index^of refraction„ 
f = 2 GHz 
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Since the ionosphere is a dispersive medium, a distinction must be made between 
group and phase velocities and group and phase indices of refraction. The lower 
atmosphere, however, is non-dispersive . Therefore phase and group velocities 
are the same and phase and group indices of refraction are the same. 


The curves in Figure 6 were obtained by ray tracing through the atmosphere to 
both satellites using (20) and (22) and the geometry of Figure 4. The range error, 
equation (5), was obtained using 


Rp - Ri + Rp +2 


C 


n dr 


V 1 - (nprp/nr)2 


(23) 
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Figure 7. Ionospheric Range Rate Error for an Actual Ephemeris 
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where Rj^ aBd E 2 are linear segments, and the integral on the right side of (23) . 
is obtained from'(l) and Snell's law for a spherically stratified' medium' . " 

nrcos0 = Up ip (24) 

where 6 is the local elevation angle of the ray at the radius r, and n is obtained 
from (20) when ray tracing through the lower atmosphere and (22) when ray 
tracing through the ionosphere . . . ■ 

For the lower atmosphere Tq is 100 kilometers and for the ionosphere r^ is • ■ 
1000 kilometers (the height of NIMBUS in this analysis). Therefore, for ray 
tracing with the ionospheric model, R 2 is zero. 


The bending angle r is given by [21] 


n. = - f 


^ dn 

cot0 — 
n 


or using (!)■ and (24) 


V = -2 npip • 10*^ 


’"o (dN/dr) dr 


/ - 

n 


2r>/l - (nprp/nr)2 


where n and N are obtained from either (20) or (22) 


The AR curves in Figure 6 were obtained from the relation 


• dAR dr 

AR = • — • — 

dr d<|> dt 


(25) 


(26) 


(27) 


where the rates d A E/dr and dr/d^ were calculated in a ray trace computer 
prograrh and the angle 0 is the total central angle subtended (Fig. 4). ‘ The rate 
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d^^>/dt was taken to be 0. 001 radians /second which is the angular velocity of the 
NIMBUS satellite-at an altitude of 1000 kilometers. This value represents an 
Upper bound for the rate of change of the central angle (in the plane of both 
satellites and the center of the earth) given the geometric configuration of the 
ATS and NIMBUS orbits. The central angle will vary with time and be less than 
this value. 

It can be seen from Figure 6 that the effects of the ionosphere are negligible in 
comparison with the effects of the lower atmosphere. Lower atmospheric bending 
at the surface is on the order of 1. 5 degrees compared with ionospheric bending 
of about 0.5 millidegrees. It is mainly the bending angle r which contributes 
to the large lower atmospheric range error of approximately 1, 7 kilometers. 
Compare this with an ionospheric range error at the surface of about -6 meters. 

The primary effect of interest in this analysis is the range rate error which is 
on the order of 100 meters/second for the lower atmosphere, hut only a few 
centimeters/second for the ionosphere. This latter follows from the fact that 
the ionospheric range error is nearly constant over the range of interest. These 
ionospheric effects are for an electron content of 10^^ electr ons /meter Since 
this quantity typically varies between 10^^ and 10^^ electrons/meter^ , the range 
rate effect should vary between about 0 . S centimeters/second and 30 centimeters/ 
second. 
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In order to further investigate the effect of the ionosphere upon' range rate, an- . 
actual ephemeris of the ATS-6 and NIMBUS-6- satellites was -used (IS, April 1976) 
along with the same ionospheric density model and lower atmospheric exponential 
refractiviiy model. , however, was. taken to be 0.42485 10^'^ electrons/ . 

meter 3 (corresponding to a total electron- content -of about 1.4 • 10 electrons/ 
meter2). Rays linking the two satellites were traced with and -without the -ionos- 
pheric model included, and the range rate along the ray path obtained by dotting 
the respective velocity vectors onto unit vectors along the ray path at each 
satellite and summing the result. The difference between the range rates as a 
function of the height of closest point of approach of the ray was then obtained. 

The results can be seen in Figure 7. The maximum effect is on the order of 
about 6 centimeters /second. 

It can be concluded from the above that ionospheric effects upon range rate are 
negligible for the frequencies involved here and for rays which pass close to the 
earth. 

At higher altitudes a simple model for ionospheric density, such as the Chapman 
model indicated in (22) should suffice for removing these effects or they could be 
removed by employing a dual-frequency system. 

If the purpose of our investigation were to recover ionospheric densities, then we 
would only be interested in rays passing through the ionosphere and not the lower 
atmosphere . In this ease the situation would be reversed and ionospheric effects 
would be much greater than lower atmospheric effects [19]. 



8.0 RECOVERY OF PRESSURE AND TEMPERATURE FROM REFRACTIVITY 
In this section it is assumed that the effects of the ionosphere and water vapor in 
the atmosphere have been removed from the recovered refractivity (either by re- 
gression techniques using some assumed model or by independent measurement) . 


Thus, all that remains is the dry refractivity N(j , a quantity directly proportional 
to density. From (2) 


Nd = 


77.6 P 
T 


(28) 


where P is pressure in millibars and T is temperature in degrees Kelvin. 


Two methods are described for calculating temperature and pressure parameters 
from dry refractivity. In both methods it is assumed that hydrostatic equilibrium 
and the perfect gas, law hold, hi the first method the dry refractivity is numer- 
ically integrated to obtain pressure directly [22] . Temperature is then deter- 
mined from Nd and P using (28). The second method involves solving for a base 
temperature T^ and temperature lapse rate j3 in a given atmospheric layer using 
values of dry refractivity in a least squares sense (weighted apriori information 
can be included). Pressure can then be calculated from the dry refractivity at 
the base of the atmospheric layer and the temperature profile. In this method 
it is assumed that temperature is a linear function of height. 


Pressure by Integration of Dry Refractivity 
The hydrostatic equation is [23] 


- * T T) A rtF. IS 


dP - -pg(h)dh (29) 

where p is the density of dry air and g(h), the acceleration of gravity, is a 
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fimctiot'of the geometric Tiei’ght h* ( [24], pg. 217). ' 


The perfect gas law is 


■ P = 


MP 

•RT„ 


. 430 )-: 


where M = the apparent molecular weight of dry air = 28.966 ([24],.p. 289) 

R = universal gas constant = 8, 31436 joules (°Kmole) ^ ( [24] , p. ,289)* 
P = pressure in millibars 
Tv -= virtual temperature in “Kelvin 
and [23] 


T = 

Ay 


T 


(31) 


' 1 - C.379- (|-) 

where T is teinperature iii degrees Kelvin, and e is the partial pressure of the i 
water vapor in the air '(in millibars) given by 


e = 


100 


7.5 (T - 273.15) 
(6.1 1) 10237-3 273.15) 


.(32) 


and'Rh.is the relative* humidity in percent. 


*The expression for g(h) is 

20^0 

g(h) = (meters/second^) 

(Jo h)2 - ' ■ 

where the local acceleration of gravity ([24], p. 488) and the effective earth’s radius ([24] , 
p. 218) are functions of latitude 

= 9.780356 [1 +0.0052885 sin20L -.(5.9)10-6 sin^ (20 l)1 

‘ ' ' 2go 

° , [(3 .08S462)'x 10-6 + (2i27) x IQ-^ cos (20J - 2 x lQ-12 cbsf40L)] 
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It is sufficient to consider g(h) constant 


g(li) = G = 9.8 meters/second^ 


(33) 


Substituting (30) into (29), and using (28) and (31) gives 

/ MG 


dP = - 


J7.6R/ VX 


—] Nddh 


(34) 


Since we are assuming no water vapor, e = 0 and (T/Ty) = 1, and (34) becomes 


dP = -0.44Nd dh 


(35) 


Integrating (35) gives 



Temperature T can now be obtained using (28) since and P are known. 


Base Temperature and Lapse Bate from Dry Refraetivity 
Equation (34) can be rewritten (again letting (T/Tv) =1) 


dP ^ 
P 



dh 


(37) 


Letting the temperature be linear with height, pressure can be obtained in closed 
form by integrating (37). 


P = Pb 


MG 



(38) 
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where and are the temperature and' pressure .respectively:-at'>tKe:base • i 
h =_hb of an atmospheric layer extending to some height h^ and j3 is the tem- 
perature lapse rate in that region. T is therefore given by 

T = Tb •+ Pfh-hb)' "foV''(hb'^h<‘hj.V (39) 


Using (38) in (28) an expression for dry refraetivity results. 




(40) 


where N 5 is the dry refraetivity at the base h = hb. 


If temperature is constaht.with height between hb -and h-j. , the corresponding 
expressions for P and Nj are 


= Pb exp 


, (h-hh) 

L RTb ^ 



and 


Nd = Ni 


MG, - , , 


for (^^b ^h^hj) 


(42) 


From (41) and (42) it is seen that a constant temperature with height gives 
f rise to';an'expohential.distribution of -pressure and'ah exponential refraetivity • 
profile. 
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Combining (40) and (42) results in an expression in two variables and ^ 
for the ratio (N^/Nt,) 


Nb 



jS 

Tb 


exp 


MG 

^b 


(h~\) 



P-0 


( 43 ) 


If three or more values of dry refractivity are available (two or more ratios), 
then least squares can be used to estimate Tj^ and j3 , 

Apriori values for the base temperature Tb and the lapse rate j3 with weighting 
can be included in the least squares technique in the following way (these apriori 
values can be obtained from a standard atmosphere or previous measurements 
such as radiosonde data or satellite radiometer data). 

Let Wgj represent the weights of the apriori values at the heights hj and 
the weights of the corresponding recovered values for the ratios of the dry re- 
fractivity (Nj /Ng). Then the following sum of squares Q is to be minimized 

with respect to and P , 

■ m 

Q = ^ wii [Fi(fo, P)-Fi(To, + w2j [(Ni/No-F,(To, P)]^ (44) 

i = 1 
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where 


Fi(To,|3) = 


[i O’! -vl '(w*‘) e*o 

-r MG' ’1 • .”■ ■ • 

exp -— (hj-ho)^ ■ iS = 0 

^ L J 


In (44) To is an apriori estimate for the temperature at the height h^, while 
Tp is the "solved for" value at hQ. is obtained from 


To - + /3(ho - hb) 


The Herglotz-Wiechert algorithm yields values of the refractive index (or 
e^ivalently refract! vity) for corresponding values of the geocentric radius r.. 
The underlying assumption in the method, however, is that refractivity is 
spherically symmetric (which follows from Snell's law). In order to obtain the 
height ho in (44), the geocentric radius of the earth Pe for the closest point of 
approach of the ray' connecting the two satellites must bfe known. As a' good' first 
approximation this can be obtained from some reference' ellipsoid for the earth 
such as" Fischer's*’ model [25] 

a^'y / 1 - 

Pe = ■ , -- - (47) 

V 1 - COs2 if 

where <p' - the geocentric latitude of the closest point of approach of the ray. 


• e = V2f - f2 
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where 


(1/f) = 298.3 
3g = 6378.155 kilometers 

The height ho is determined from 

~ ^0 “ Pe 

where Tq is one of the radii obtained by the inversion technique. 


( 49 ) 


(50) 


The value obtained for Pg from (47) does not represent the geoid surface (a 
surface of equal potential) exactly. However, a correction term can be added 
to pg in (47) in order to more closely approximate the geoid surface [26] . 

It should be noted that althou'gh Tq and jS are the parameters of a linear model 
for temperature, the least squares procedure is nonlinear in these variables. 
For computation purposes, a direct approach to solving for these parameters 
is to use a search technique to minimize Q in (44), This is not costly In terms 
of computation time inasmuch as there are only two parameters. In addition 
there is no need to calculate the partial derivatives. 


Since any numerical procedure is subject to errors, both methods described 
above have certain limitations. With the numerical integration scheme a 
sufficient number of values of dry refractivity over a wide height range must 
be available to obtain pressure accurately. The main advantage, however, of 
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this method is that no assumptions are made concerning the temperature struc- 
ture. On the other hand, if there are enough points in a given height range and 
an assumption of a linear temperature structure is not unreasonable, then least 
squares can be used to obtain the temperature profile. Instead of just one layer, 
the lower atmosphere could be broken up into two or more regions with differing 
lapse rates. This would, of course, necessitate more data points (values of 
dry ref r activity). 

9.0 ERROR ANALYSIS FOR RECOVERED PARAMETERS 
In this section an error analysis is performed in order to assess the sensitivity 
of random and bias type errors on the recovered parameters. The error sources 
considered are measured range rate, position and velocity of the satellites (earth- 
fixed coordinate system) , -and height. For simplicity and ease of computation, an 
exponential model was chosen as the refractivity model for the. lower atmosphere. 
It was felt. that .the errors from a more sophisticated model would be of the same 
order of magnitude for the same size error sources. .The exponential model is 
given by • 

N(j == exp(ai’ + E2 s) (51) 

where s = (0, 01)h.- 1, h =• height (kilometers), and ai = -9.42, &2 = -15.22km. 

The nominal ephemeris used was an actual ephemeris of the ATS-6 and NIMBTJS- 
6 satellites on 15 April 1976 [27] . 
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Using Equations (1), (4), (24), and the model in (51), rays were traced between 
the positions of the satellites (as determined from the ephemeris) at five second 
intervals for a total of six data points. The range rate along the ray path was 
then calculated (from velocity components of the satellites) and served as a 
nominal "observable. " The position and velocity components of the satellites 
were perturbed (either randomly or by adding a bias) and the parameters 
and solved for in the least square sense using the "observed" range rate and 
the perturbed ephemeris and employing ray trace procedures. In the case of 
simulated errors in the range rate, the perturbed range rate became the "ob- 
servable" and the nominal ephemeris was used. Then, as before, and'a 2 
were solved for using least squares and ray tracing. 

A Monte Carlo approach was used to simulate the effects of random errors in 
position, velocity, and range rate. A random sample of ten sets (each set con- 
taining six data points spaced five seconds apart) was used for the simulation 
and statistics taken on this sample. The standard deviations of the errors in 
position and velocity of each satellite were taken as being representative of 
ATS-6 and NIMBUS-6 short term accuracy [17]. hi the case of position errors 
both satellite position components were perturbed and in the case of velocity 
errors both satellite velocity components were perturbed. 

For any perturbed set of conditions, after solving for the parameters and ^^2 
in the least squares sense, there result errors Aaj and A &2 from the nominal 
conditions and a. 2 . 
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^^2 ~ ^2 ~ ^2 


( 52 )-- 


Neglecting for the moment errors in the height h, an error in N^ can be approxi- 
mated by 


or 


ZaNA /3N<3 

AN(i = I — ) Aai + 

[daj [da^ 


Aai 


ANd 

Nd 


Aa^ "J" sAa2 


(53). 


(54) 


for bias type errors. 


In the case of random errors, Aa^ and Aa 2 are random variables and.the ya;ri- , 
ance of (54)-becomes 

^ = .(56) 

where cr^a^> °^Aa 2 * ^Aaj^Aa 2 variances of Aaj^ and Aa 2 and the 

covariance between Aa^ and Aa 2 which can be estimated from the sample. 


Using (36)_ and (51J the pressure can be written as 

-44 

P = — exp ('a. +■ ajs) (56) 

^2 

For pressure the following bias and random error expressions result. 
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s 


J_ 

^2 


(57) 


P 


= Aaj 



Aa2 


and 



+ 2sa 


Aaj^ Aa2 


_t_ q2 ^2 

^ 5 ^Aa, 


(58) 


For temperature, using (28) and (51) 


and therefore 


and 


3414.4 

^2 



(59) 


(60) 


(61) 


for bias and random type errors respectively. 


For errors in height, using (51) we can write 


and 


ANd 

Nd 


(0.01) aj Ah 



= 10 


r4 


^2 


(62) 


(63) 


Since N is exponential, temperature is constant, and therefore there is no 
error with height. 
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Results of the error analysis can'be seen in Figures 8 through 13 . It is apparent 

i 

from the figures that the percent error in recovered refractivity and pressure is 
linear with parameter error for a given height. Also, recovery of the param- 
eters is more sensitive to errors in velocity of the satellites than to position 
errors, in particular, the position errors in the geostationary satellite. 

It should be noted that the geoid varies approximately -110 to +80 meters over 
the globe with a standard error of about 5 meters from a reference ellipsoid 
[26] , Thus, using Equation (62) with an error of 5 meters and the value of'a 2 
as^given in (51), recovered refractivity is in error less than 0.1%. 

10.0 AN EXAMPLE 

Figure 5 shows the lower atmospheric effect upon, one-way range rate between 
the NIMBUS-6 satellite and the ATS-6 satellite on 15 April 1976, Radio occulta- 
tion of the tracking signal between the two spacecraft occurred at about 20 hours; 
47 minutes, 0 seconds GMT, and lasted for approximately half a minute, NIMBUS 
(in a near-circular almost polar orbit of 1100 kilometers height) was on a northern 
trajectory at about 64“ South latitude and 116“ West longitude while ATS was geo- 
stationary (equatorial plane) at 35“ East longitude. The closest point of approach 
of‘the ray connecting the satellites was at 74° South latitude and 20“ West longi- 
tude, The nearest radiosonde station at the approximate latitude was McMurdo 
Sound station^located at 77“ 53”* South and 116“ 44™ -East and the radiosonde data-. 
was taken at OGMT on 15 April 1976. 
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NIMBUS ACCURACY 


X - 15 meters x - 5 cm/sec 

Y - 10 meters y - 3 cm /sec 

z - 20 meters z - 8 cm/sec 

ATS ACCURACY 

X - 50 meters x - 0.2 cm/sec 

y - 50 meters y - 0.2 cm/sec 

2-50 meters z - 0.5 cm/sec 

RANGE RATE ACCURACY 
0.5 cm/sec 


RANGE 



Figure 8. Effect of Random Errors on Recovered Refractivity (Monte Carlo) 
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NOTES: 

- • EXPONENTIAL REFRACTIVITY USED AS 
NOMINAL FOR RAY TRACING 

ATS-6/NIMBUS-6 EPHEMERIS OF 15 
APRIL 1976 used" AS'NOMINAL 

' • % ERROR IN REFRAcfiyiTY IS 

LINEAR WITH PARAMETER ERROR 
FOR A GIVEN HEIGHT . 



Pignre 9. Effect of Bias lype Errors on Recovered Refractivity 




NIMBUS ACCURACY 

X - 15 meters x - 5 cm/sec 

y - 10 meters y - 3 cm/sec 

z - 20 meters z - 8 cm/sec 

ATS ACCURACY 

X - 50 meters k - 0.2 cm/sec 

y - 50 meters y - 0.2 cm/sec 

z - 50 meters 2 - 0.5 cm/sec 

RANGE RATE ACCURACY 
0.5 cm/sec 


RANGE 



ERROR IN RECOVERED PRESSURE (%) 

Figure 10. Effect of Random Errors on Recovered Pressure (Monte Carlo) 
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NOTES: 

• EXPONENTIAL REFRACTIVITY USED AS 
NOMINAL FOR RAY TRACING 

• ATS-6/NIMBUS-6 EPHEMERIS OF 15 
APRIL 1976 USED AS NOMINAL 

• % ERROR m PRESSURE IS 
LINEAR WITH PARAMETER ERROR 
FOR A GIVEN HEIGHT 



Figure 11. Effect of Bias Type Errors on Recovered Pressure 



ERROR TYPE 

% ERROR IN TEMPERATURE 

Position 

0.20 

Velocity 

0.68 

; Range Rate 

1 

0.05 


NIMBUS ACCURACY 


X 

- 15 meters 

• 

X - 

5 cm/sec 

y 

- 10 meters 

• 

y - 

3 cm/sec 

z 

- 20 meters 

9 

Z - 

8 cm/sec 


ATS ACCURACY 


X 

- 50 meters 

• 

X - 

0.2 cm/sec 

y 

- 50 meters 

« 

y - 

0.2 cm/sec 

X 

- 50 meters 

* 

z - 

0.5 cm/sec 


NOTES: 

• Exponential refractivity used as nominal for ray tracing and therefore the nominal temperature 
distribution’is a constant for all heights. 

• ATS- 6/N I MB US-6 ephemeris of 15 April 1976 used as nominal. 

Figure 12, Effect of Random lYpe Errors on Recovered 
Temperature (Monte Carlo) 
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ERROR SOURCE 

% ERROR IN TEMPERATURE 

NIMBUS POSITION 

-0.08 

Ax = Ay = Az = 20 meters 


ATS POSTIION 

0.03 

Ax = Ay_ - Az. = 10-kilorneters 

. 

NIMBUS VELOCITY 

-0.14 

Ax = Ay = Az = 1 cm’/second 


ATS VELOCITY 

0.13 

Ax = Ay = Az = 1 cm/second 


RANGE RATE 

-0,009 

Ar = 1 mm/sec 



NOTES: 

• EXPONENTIAL REFRACTIVITY USED AS NOMINAL FOR RAY TRACING 
AND THEREFORE THE.NOMINAL TEMPERATURE DISTRIBUTION IS A 
CONSTANT FOR ALL HEIGHTS 

• ATS-6/NIMBUS-6' EPHEMERIS OF 15 APRIL 1976 USED AS NOMINAL 

• % ERROR IN RECOVERED TEMPERATURE IS LINEAR WITH 
PARAMETER ERRORS 

Figure 13. Effect’ of Bias Type Errors on Recovered Temperature 


The ephemeris of both satellites during occultation was obtained by propagating 
their state vectors at epoch [27-] , (a time .point close to occultation) through this 
period of time. The states at epoch were determined from non-occulted tracking 
data. In the case of ATS-S approximately 3300 observations of two-way range 
and range rate tracking measurements made from the Madrid tracking station, 
at a data rate of one per ten seconds were used in the solution. For NIMBUS-6 
about 2500 measurements of the 4-way range and range rate (see Section 4. 0)' 
from Madrid were used in the solution. An 8x8 gravity field for NIMBUS and 
a 4 X 4 gravity -field! for ATS were used also. Drag was not applied for either 
satellite. However, solar pressure was used in the solution for ATS-6 [273. 
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Refractivity was recovered using the ephemeris of the satellites and all but the 


last nine points of data in Figure 5 by the Herglotz-Wiechert direct integral 
transform technique and by a least squares procedure as a check. The last 
nine points of data (from 20 hours, 47 minutes, 25 seconds, to 20 hours, 47 
minutes, 33 seconds) were excluded from the calculations as it was felt that 
these points were quite noisy, possibly due to loss of lock of the tracking signal 
or multipath effects or a combination of both.* Results of the* recovered re- 
fractivity and refractivity calculated from radiosonde data are shown in Figure 
14. Recovery by both methods differed by a maximum of about 2. 3% at 10. 5 
kilometers. 

Retrieval of pressure and temperature from dry refractivlly was obtained by 
the two methods described in Section 8.0. Figure 15 shows pressure calcula- 
tions against radiosonde. Direct integration of dry refractivity agreed more 
closely with radiosonde data with a maximum of about 3. 5% at 10. 5 kilometers. 
Temperature recovery is compared with radiosonde data in Figure 16. The 
solid line shows temperature for a standard atmosphere. Also shown are re- 
covered temperature profiles using least squares with and without apriori 
weighting (the standard atmosphere was used as an apriori profile) and tempera- 
ture obtained from dry refractivity and pressure (via numerical integration of 


*When these points were included in the solution the resulting radii of closest approach of the ray 
connecting the two satellites did not decrease monotonically with time (i.e,, as NIMBUS went below 
the horizon). This should occur for Snell’s law to hold. 
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ATS/NIMBUS OCCULTATION 15 APR.1L 1976 



Figure 14, Comparison of Refractivily Recovery with 
Refractivity Calculated from Radiosonde Data 



Fi^re l5. Pressure Recovery Versus Radiosonde Data ’ 

1 . 
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ATS/NIMBUS OCCULT ATIOISl 15 APRIL 1976 



TEMPERATURE FROM INTEGRATED 
DRY REFRACTIVITY AND DRY 
REFRACTIVITY 

APRIORI PROFILE {STANDARD 
ATMOSPHERE) 

OCCULT ATION DATA WITH WEIGHT 
OF 1 AND APRIORI PROFILE WITH 
A WEIGHT OF 1 

OCCULTATION DATA ONLY 


Figure 16. Temperature Recovery Versus Radiosonde Data 


dry refractlvity). Without weighting the maximum deviation of the least squares 
profile was 1.5% at 10. 5 kilometers. 


The radius of the geoid at the latitude and longitude of the point of closest ap- 


proach of the ray was subtracted from the calculated radii of closest approach 
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using Fischer's ng.p4el and. a correction term for the .geoid [26l to obtain height. 
Since .the recovered data was at an altitude above 10 kilometers water vapor 
effects were not included. The local time of the pass was about 2 a,m, and 
therefore ionospheric effects were considered minimal and were not included, 

11,0 DISCUSSION AND CONCLUSIONS 

This analysis has demonstrated the feasibilily of recovering refractivity profiles 
and pressure and temperature distributions for the lower atmosphere from radio 
occultation data taken during satellite- to-satellite tracking at S-band frequencies. 
One satellite, ATS-6, is in geostationary orbit (equatorial plane) and the other 
is in a near circular polar orbit. The only assumption made in the inversion 
technique is that refractivity is spherically symmetric (Snell's law holds). 

Two methods have been investigated for determining pressure and temperature 

i 

4 V 

distributions from the recovered refractivity., In both methods it is assumed 
that the hydrostatic equation and perfect gas law hold. In a dry atmosphere re- 
fractivity is directly proportional to density, a function of both temperature and 
pressure. Using this fad, the first method consists of numerically integrating 
dry refractivity to obtain the pressure directly, without any assumptions con- 
cerning the temperature structure. Temperature can then be determined from 
pressure and dry refractivity. In the second method least squares is used to 
obtain the parameters of a linear temperature profile with height over a given 

2 • H 

region of the lower' atmosphere (temperature at the base of the level and the 
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lapse rate in that region). Apriori information with appropriate weighting (such 
as a standard atmospheric temperature profile or previous radiosonde or satel- 
lite radiometer measurements) can be included. Pressure is then obtained from 
the recovered temperature profile and dry refractivity at the base of the layer. 
Both methods yield identical results for a linear temperature profile in a dry 
atmosphere. 

At low altitudes the effects of water vapor in the atmosphere must be taken into 
account, either by independent measurement or by regression techniques using 
some assumed model [18]. At these altitudes and particularly at S-band fre- 
quencies (the carrier frequency of the tracking signal in this analysis is 2 GHz) 
ionospheric effects are almost negligible compared to lower atmospheric effects. 
At higher altitudes water vapor does not present a problem, but ionospheric 
effects increase somewhat. These effects can be removed by employing a dual 
frequency system. In the case of a single frequency system a simplified ionos- 
pheric model for electron density should suffice to remove these effects. 

Results of an error analysis, where a simple exponential refractivity model for 
the lower atmosphere is assumed, indicate that recovery of atmospheric param- 
eters is more sensitive to errors in the velocity components of the satellites than 
to errors in the position components of the satellites. This is particularly true 
for the geostationary satellite where a bias -type error as great as 10 kilometers 
in each of the position components gives rise to an error of less than 0.2% in 
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recovered refractivify up to approximately 20 kilometers. Errors in the velocity 
components of each satellite on the order of 1 centimeter/second give rise to 
errors less than 0, 8% in pressure and refractivity and about 0. 1 % in temperature 
up to about 20 kilometers. 

An example has been given demonstrating refractivily recovery and calculation 
of pressure and temperature distributions from radio occultation data obtained 
during satellite-to-satellite tracking of the NIMBUS-6 satellite (in near circular 
polar orbit) by the geostationary satellite ATS-6, An exponential refractivity 
profile was obtained by least squares procedures and used as a check against 
recovery by the Herglotz-Wiechert inversion algorithm. Recovery by both 
methods differed by a maximum of about 2. 3% at a height of 16 kilometers. 

Both methods for calculating pressure and temperature were used and the re- 
sults compared against radiosonde data. In the case of pressure recovery, 
direct integration of dry refractivity agreed more closely with radiosonde data 
with a maximum deviation of about 3,5% at 10.5 kilometers. The least squares 
procedure with and without apriori information yielded a temperature profile 
closer to the radiosonde data with a maximum deviation of about 1. 5% at 10. 5 
kilometers. 

With the geometrical configuration of the satellites as analyzed here, different 
levels of the atmosphere can be probed (as opposed to two satellites in the same 
orbit with a fixed spacing between them). 



For a fixed inclination occultations will occur a few degrees in latitude below the 
value of inclination both in the northern and southern hemisphere, ha addition the 
occultations will occur at different longitudes. For example, a satellite In a 
1000 kilometer orbit about the earth will make about 14 revolutions per day which 
means 28 occultations at approximately the same latitude but varying longitudes. 
Thus, greater coverage could be accomplished by having a number of satellites 
at different inclinations. No additional electronic hardware (other than that re- 
quired for satellite-to-satellite tracking) would be required onboard the space- 
craft, The only requirement would he for the low satellite to be tracked as it 
went down or came up over the horizon with respect to the geostationary satel- 
lite at a sufficiently high data rate to ensure enough points for inversion pro- 
cedures. Such a system could be used as an adjunct method for monitoring 
world-wide weather patterns . 
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APPENDIX A 

HEEGLOTZ-WIECHERT INVERSION* 


Referring to Figure 4, let r^ be the radius of the outer limit of the atmos- 
phere. Then, the electrical range from N (ISflMBUS) to A (ATS) can be broken 
into three segments: the straight line distance outside the atmosphere from N 
to D, the electrical range along the ray path (inside the atmosphere) from D to 
B, and the straight line distance outside the atmosphere from B to A. The 
range from D to B is 


R = 2 


f "£) 


(A.1) 


where n = n (r) is the index of the refraction at radius r, i is the angle of in- 
cidence, ds is an incremental path length, and rp is the radius of closest ap- 
proach of the ray to the Earth. 


Using Snell's law for a spherically stratified medium, a ray parameter p 
called the impact parameter is defined 

p = n.r sin i = sin ip (A.2) 

The parameter p is a constant along the ray path, and, in the absence of an 
atmosphere is the closest point of approach of the ray to the Earth. 


• *This derivation follows that in [11]. 
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Defining a. new variable rj , 


T? = n r = 


sin 1 


Equation (A. 1) can be written 


! V^\ / dr\ 

R = R (p) = 2 I (— ) ( — ]dT? 

Jp W / \ 


d-qj 


Set 


2 2 
y = - to 


1 1 
W = P - Tq 


Then (A. 4) becomes 


/ d In r \ 

R(w) ^ J dy 


In order to solve (A. 6) for In r the Abel transform pair is used 

J rw 

I s'(y) (y - dy 

o 


■I 


g(y) = (kvr)"^ I f(w) (w - y)"’'^" dw 


Equation (A< 6) is in the form of the first equation in (A. 7) with f = F 
k = 2, and 


g'(y) = (y + rj) 


d In r 
dy 


A~2 


CA.3) 


(A- 4 ) 


(A. 5 ) 


(A. 6) 


(A.7) 


(A-8) 
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Differentiating the second equation in (A. 7) and changing back to the vari- 
ables 7} and p 


V 


ir, I 


R(p)(p2-T?^r’'^(2p)dp)- (A.9) 


di? 


Integrating (A. 9) from ?j to 


In 




Integrating (A. 10) by parts gives 




,2 


■e/:-f 


dp 


pR(p) (p2 dp du 


(A.ll) 


Interchanging the order of integration in the second term of (A. 11) gives 


0/:-f 


PR(P)(P^ dpdu 




pR(p) dp I u ^ (p^ -u)“^ dn 

'17 


(A. 12) 


0 f 


Vp^ - 77^ 


pR(p)i ' ; + (— 1 hi 

2p^77^ \2p- 


(i) 


P+Vp^-1?^ 


►dp 


A-3 
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Using 

cosh“^ (z) = J* (t^ - 1)"''^" dt = In [z + (z^ - I)’''" ] 
tBe second term of (A. 11) becomes 

pR(p) (p^ - u^)"^ dp du 


m-'-r 


•r. 



PR(P) 


V 


1 /y\ , / P ^ 

— (P^ ■^(t) cosh'M - 

,p2 \p3-/ \7jyj 


dp 


(A. 13) 


(A. 14) 


Equation, (A. 10) can now be written 

/r(r?)> 

In 



(p2 - - f—] (p^ - T?^ 


dp 


or 


In 




(A.15) 


Using 


A /n cosh- 

dp [Vp/ \n 


;)} ■ -£) 


cosh”^ /-T^— = 

P"/ \vj \p/\VP -• 


and partial integration, Equation (A. 15) is 


/r(T?)\ 1 

In( } = — I cosh"^ I — )( — Jl — Idp 

\ Tq / 1 T \v/\p/\dp/ 

O 


(A- 16) 


(A-17) 
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At the closest point of approach of the ray r 

r 


r 


p 




(A. 18) 


Using (A. 2) 


/r^smiA fl /pwlWdRN 'I 

n =[ ^ 1 = (siniQ)expJ— f cosli"^ (— )( — 1{ — Idpl 

\ rp / [^Ji . V7?Ap/Vdpy J 


(A.19) 


The range along the ray path from D to B, R , can be written as the sum of 
the range in the absence of an atmosphere, R^ , plus the range residual due to 
an atmosphere, 

R=Ry+Rg • (A.20) 


In the absence of an atmosphere (A. 19) becomes 


1 = (sin ig) exp 




(A.21) 


Substituting the expression for sin i^ obtained from (A.21) into (A.19) and 
using (A.20) gives 


np(p) = exp 





(A.22) 


In the absence of an atmosphere 

Ru = 2pcoti^ 


(A.23) 
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Therefore 


R, = R-R„ 



n-d s - 2p cot Iq 


(A.24) 


Using the electrical range from N to A and the fact that the range from N 
to D and the range from B to A are straight lines can be written 

= Re - VRi-p2 - (A.25) 

where Rj^j and R^ are the radius vectors of the low and high satellites 
respectively. 

Then, differentiating (A.25) gives 
dRj, . / dt \ 

— = Re ( — ) + p[(Ri -p2)-^^ + (R|, -p2)-/n (A.26) 

dp \ dp / 


The expression in (A. 26) is used in (A. 22) and contains the measured one 
way. range rate along the ray path from NIMBUS to ATS. 
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